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ABSTRACT 

We present a model of optically thin, two-temperature, accretion flows using an exact 
Monte Carlo treatment of global Comptonization, with seed photons from synchrotron and 
bremsstrahlung emission, as well as with a fully general relativistic description of both the 
radiative and hydrodynamic processes. We consider accretion rates for which the luminosities 
of the flows are between ~ 10"^ and 10"^ of the Eddington luminosity. The black hole spin 
parameter strongly affects the flow structure within the innermost ^10 gravitational radii. The 
resulting large difference between the Coulomb heating in models with a non-rotating and a 
rapidly rotating black hole is, however, outweighed by a strong contribution of compression 
work, much less dependent on spin. The consequent reduction of efi'ects related to the value 
of the black spin is more significant at smaller accretion rates. For a non-rotating black hole, 
the compressive heating of electrons dominates over their Coulomb heating, and results in 
an approximately constant radiative efficiency of a 0.4 per cent in the considered range of 
luminosities. For a rapidly rotating black hole, the Coulomb heating dominates, the radiative 
efficiency is ^ 1 per cent and it slightly increases (but less significantly than estimated in some 
previous works) with increasing accretion rate. Our study neglects the direct heating of elec- 
trons, which effect can lead to larger differences between the radiative properties of models 
with a non-rotating and a rapidly rotating black hole than estimated here. Flows with the con- 
sidered parameters produce rather hard spectra, with the photon spectral index F ~ 1.6, and 
with high energy cut-offs at several hundred keV. We find an agreement between our model, 
in which the synchrotron emission is the main source of seed photons, and observations of 
black-hole binaries in their hard states and AGNs at low luminosities. In particular, our model 
predicts a hardening of the X-ray spectrum with increasing luminosity, as indeed observed 
below ~ O.OILe or so in both black-hole binaries and AGNs. Also, our model approximately 
reproduces ffie luminosity and the slope of the X-ray emission in Cen A. 
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1 INTRODUCTION 

Optically thin, two-temperature accretion flows have been consid- 
ered as an explanation of a variety of black hole systems and a 
substantial work has been done for investigation of their dynami- 
cal and spectral properties (see, e.g., Narayan & McClintock 2008 
for a review). However, the developed models still involve several 
approximations which significantly reduce their accuracy. These 
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involve the use of a pseudo-Newtonian potential of Paczyriski & 
Wiita (1980), which fails in the innermost region (particularly if 
rotation of the black hole is considered), where most of the grav- 
itational energy is dissipated, as well as local approximations of 
Comptonization, which appear to be particularly incorrect in opti- 
cally thin flows (see Xie et al. 2010; hereafter XIO). Previous at- 
tempts to improve these two weaknesses are discussed in Section 

In this paper, we extend our previous treatment of global 
Comptonization from XIO by using the hydrodynamical model 
from Manmoto (2000; hereafter MOO), and develop a self- 
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consistent model involving a fully general-relativistic (GR) de- 
scription of both the hydrodynamical and the radiative processes. 
We consider moderate values of accretion rate, at which Coulomb 
coupling between ions and electrons is relatively weak compared to 
the viscous heating of ions and, therefore, the ion temperature (and 
hence the total pressure) is not affected significantly by the details 
of the description of radiative processes. This allows for a slightly 
simplified treatment of the flow structure in our computations lead- 
ing to a self-consistent solution. 

Already the foundational papers proposed that tenuous, two- 
temperature flows may be responsible for the hard spectral states of 
black-hole binaries (Ichimaru 1977; see also Narayan & Yi 1995) 
and for the low nuclear luminosities in radio galaxies with large 
radio lobes (Rees et al. 1982). We consider two ranges of the key 
parameters (black hole mass and accretion rate) which should be 
relevant for these two major application areas. For each of these two 
cases, we illustrate the impact of the black hole spin by considering 
a non-rotating and a maximally-rotating black hole. On the other 
hand, we neglect here the dependence on some other parameters of 
hot-flow models. 

In particular, we consider only weakly magnetized flows, with 
the magnetic pressure of 1/lOth of the total pressure, and we ne- 
glect the direct viscous heating of electrons; see Section 13.61 for 
a discussion of these two assumptions. Furthermore, we neglect 
outflows, which may play an important role in some systems (e.g. 
Yuan, Quataert & Narayan 2003). Then, our solutions are strongly 
dominated by the advection of energy by ions, as in the original for- 
mulation of the advection-dominated accretion flow (ADAF) model 
by, e.g., Narayan & Yi (1995). 

In this paper, we focus on modelling an innermost part of an 
accretion flow, namely inside 10'i?g, where R^ = GM/c^ is the 
gravitational radius and M is the black-hole mass, where the bulk 
of the observed radiation is produced. Global Comptonization has a 
net effect of Compton heating outside ~ lO'^R^ (e.g. Park & Ostriker 
2001, Yuan, Xie & Ostriker 2009, Yuan & Li 201 1), which efl'ect is 
not considered here. 



2 THE MODEL 

We consider a black hole, characterised by its mass, M, and angular 
momentum, J, surrounded by a geometrically thick accretion flow 
with an accretion rate, M. We define the following dimensionless 
parameters: r = R/Rg, a = JI{cR„M), m = M/Me, where Me = 
Le/c^ and Le = 4nGMmpc/crT is the Eddington luminosity. The 
inclination angle of the line of sight to the symmetry axis is given 

by/^obs = COS^obs- 

We refine here our previous study presented in XIO by in- 
cluding a fully GR hydrodynamical description of the flow. Our 
GR hydrodynamical model follows strictly MOO, except for a small 
difference in the advective term, as noted below equation Q. We 
have thoroughly tested our hydrodynamical model by comparing it 
with the numerical code and results of Li et al. (2009), which work 
also follows MOO. Other differences with respect to XIO involve: 
(i) the vertical structure of the flow: we assume here that the tem- 
perature and velocity components are vertically uniform, while the 
density distribution is given by p(R,z) = p(R,0)exp(-z^/2H^), 
where H is the scale height at r; in XIO we assumed a tangentially 
uniform structure of the flow; 

(ii) the direct viscous heating of electrons, which is neglected here 
and included in XIO; 



(iii) an outflow, neglected here and included (with a large magni- 
tude) in XIO. 

We consider the pairs of the main parameters of {M = 2 x 
10** M0, m = 0.1) and (M = IOMq, m = 0.5), which correspond 
to low-luminosity AGNs and hard states of black-hole binaries, re- 
spectively. For each pair of (M, m) we consider two values of the 
spin parameter, a = and a = 0.998. We assume the viscosity pa- 
rameter of Q- = 0.3 and the ratio of the gas pressure (electron and 
ion) to the total pressure of ySs = 0.9. The latter parameter deter- 
mines the strength of the magnetic field in the accretion flow, and 
the synchrotron emission can then be determined. For it, we follow 
the method of Narayan & Yi (1995), which was applied also by 
MOO and XIO. 

We define the compressive heating rate of electrons per unit 
area (in the local reference frame) of the flow, 



A„ 



Mpe dlnp 
'IttRp dR ' 



(1) 



and the rate (per unit area) of advection of the internal energy of 
electrons. 



ein 



Mp, 



dlnr. 



2nRp(r^ - 1) dR 



(2) 



where p is the mass density and p^ is the electron pressure in the 
midplane, Fe is the effective adiabatic index [see, e.g., eq. (62) 
in MOO for the definition]; for simplicity of notation we hereafter 
skip the 'e' subscripts in symbols denoting electron heating/cooling 
rates. Note that in our definition M > 0. Then, A^ompi is always 
positive (therefore, we define it as a heating term) as the density 
gradient is always negative in our models. In all models, Q]„f is 
negative (representing the release of internal energy) in the inner- 
most region, at larger distances it is positive (the internal energy 
is stored) or vanishes. Following the usual notation, we define also 
the electron advection rate, Qadv - Qmt - A-compi. which is typically 
negative (i.e. the advective effects result in effective heating) ex- 
cept for an outer (beyond r^, see below) region in the local models. 
The above forms of Acompi and Qint follow directly from eqs. (54) 
and (60) in MOO. We do not use the simplifying approximation, 
dln/f/dlni? = 1, adopted in MOO; we checked that it yields devi- 
ations of the advective term from the actual values by a factor of a 
few. 

With the above definitions, the electron energy equation is 



— Aie + At;omp|- — Qrdd ~ Qirtt^ 



(3) 



where Qrad is the radiative cooling rate per unit area (including syn- 
chrotron, 2synch> and bremsstrahlung, Qbrcm. and their Comptoniza- 
tion, 2compt) and Aie is the electron heating rate per unit area by 
ions via Coulomb collisions. 

The values of Aje, Acompr and Qr^i integrated over the whole 
body of the flow are denoted by Aje.tot. A;.o„,p,.tot and Qiad.tot. respec- 
tively. Due to the effect of capturing by the black hole and the grav- 
itational redshift (see Section [3.3.3l for details), the total power ra- 
diated by the flow, 2rad.tot> is higher than the luminosity, L, detected 
far away from the flow. 

For completeness we give also the ion energy equation: 



= gdiss + Aco 



QinLi i 



(4) 



where the compressive heating and the advection of the internal 
energy of ions, respectively, are given by 



A __M£i_dln£ 



Mp, dlnT, 
'inRpiT^-l) dR ' 



(5) 
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and the viscous dissipation rate, per unit area, is given by 



Qdiss = - 



rjA 



dfi 



, (InY^apH—, 
r' ar 



(6) 



wiiere p = {p\ + Pc)IPb., Pi is tlie ion pressure, Fi is the ion adiabatic 
index, fi is the angular velocity of the flow, y^ is the Lorentz factor 
of the azimuthal motion and A = r* + r^a^ + 2ra. The above form 
of 2diss results from the usual assumption that the viscous stress is 
proportional to the total pressure, with the proportionality coeffi- 
cient a; cf. eqs. (48) and (54) in MOO. This is almost certainly an 
oversimplified approach; see Gammie & Popham (1998) for a dis- 
cussion of various issues related to implementing a proper descrip- 
tion of the viscous stress in GR models. The magneto-rotational 
instability (MRI) is now widely believed to be an origin of the vis- 
cous stress and numerical simulations of this process (e.g. Krolik, 
Hawley & Hirose 2005) show radial profiles of the stress which 
are different from analytic prescriptions considered in the accretion 
theory. The uncertainty regarding the form of the stress is the major 
shortcoming in attempts to develop an analytical description of ac- 
cretion flows. The form of the viscous stress, oc ap, is assumed here 
for computational simplicity. It has an unphysical property of lead- 
ing to a negative dissipation rate close to the event horizon in some 
cases. However, this particular weakness of our model appears to be 
not important for our results. Namely, we find that Q^i^^ < occurs 
at r < 3 only in our models with a = 0; in these models the heating 
of ions is strongly dominated by the compression work at r < 10. 
Therefore, details of the viscous heating of ions do not affect the 
flow structure. Clearly, the accuracy of Qdiss would be much more 
important in models taking into account direct viscous heating of 
electrons. 

Our modelling of Comptonization makes use of a Monte Carlo 
(MC) method taking into account all effects relevant for the Kerr 
metric (Niedzwiecki 2005; see also Niedzwiecki & Zdziarski 2006 
and XIO). The seed photons are generated from synchrotron and 
bremsstrahlung radial emissivities of the flow, taking into account 
the vertical profile of the density; their transfer and energy gains in 
consecutive scatterings are affected by both the special relativistic 
and gravitational effects. 

The MC Comptonization simulation and the solution of struc- 
ture equations are treated separately in our model. Therefore, we 
have to iterate between their results to find a mutually consistent 
solution. In principle, a full set of hydrodynamical equations should 
be solved after each MC simulation, as in XIO. However, in the cur- 
rent study we apply a slightly simplified procedure, motivated by 
the fact that the structure of the flow is fully determined by the ion 
energy equation (this condition is, however, satisfied only within 
inner several hundred R^, see below), while the electron energy 
equation determines the electron temperature. Under such condi- 
tions, our procedure of refining the description of Compton cooling 
should affect only the electron temperature and not the other flow 
parameters {T\,p, H, and the velocity field). Then, we first solve the 
dynamical structure with a local approximation of Compton cool- 
ing and then we iterate between the MC Comptonization results and 
the solutions of the electron energy equation, keeping the flow pa- 
rameters, except for T^, unchanged. After each MC simulation, we 
determine the radially-dependent Compton cooling. We then use 
this rate to solve the energy equation ([3]l, taking into account the de- 
pendence of all relevant processes (synchrotron and bremsstrahlung 
emissivities. Coulomb exchange rate, advective terms) on T^. We 
repeat the above procedure until it converges. 

For the initial solution, the outer boundary of the flow is set at 
r = 2 X lO'* (see section 2.4 in MOO for details). The MC simula- 



tions are performed within r„ut = 1000, using the hydrodynamical 
solution at r < rout ', note that details of this inner part of the solution 
should not depend on the outer boundary condition (see e.g. fig. 5 in 
Narayan, Kato & Honma 1997). For our final solution, we balance 
the electron energy equation up to r ^ 700. Due to boundary effects 
in MC simulation, combined with strong sensitivity of the advective 
term on the temperature gradient, the solution is very unstable at r 
close to Tout. Furthermore, we note that an assumption underlying 
our computational procedure fails in these outer regions. Namely, 
all our final solutions are characterised by an electron temperature 
larger (typically by a factor of ~ 3 at r < rout) than the initial so- 
lutions and the electron pressure can give a non-negligible contri- 
bution beyond several hundred R^. In that region, the scale height 
would increase, which effect is not taken into account here; how- 
ever, we expect it to be a rather minor, < 10 per cent, effect - such 
a change of H/R was found, in response to a similar in magnitude 
(but opposite) change of T^, in XIO. For one set of parameters, 
namely M = IOMq, m = 0.5 and a = 0, T^ approaches T^ at 
r ^ 1000, so the flow loses a two- temperature structure. We skip 
here more careful investigation of the outer region because it gives 
a negligible contribution to the observed radiation. Furthermore, a 
detailed study of the outer region should involve a proper descrip- 
tion of additional, poorly understood effects (like a mechanism of 
transition from an optically thick disc to a tenuous flow; e.g., in 
hard states such a transition is likely to occur at r < 1000) which 
are beyond the scope of this paper. 



3 RESULTS 

Figs [T] and |2] show the change of the heating and cooling rates, and 
the electron temperature, between our initial and final solutions. 
Figs [3] and |4] show the remaining parameters of the flow - these 
are assumed to be the same in the initial and final solutions. We 
can clearly see how the value of the spin parameter, a, affects the 
flow properties. Rotation of the black hole stabilizes the circular 
motion of the innermost part of the flow (an effect directly related 
with properties of test particle motion, analogous to the well-known 
dependence of the innermost stable orbit on a in Keplerian discs) 
which yields a higher density (through the continuity equation). 
This, in turn, implies both a stronger Coulomb coupling (see Figs[Tf 
and|2f) and a larger optical depth (see Figs[3^ and|4^) in the inner- 
most region in models with a = 0.998. On the other hand, heating 
by the compression work is typically stronger than, or comparable 
to, heating by the Coulomb energy transfer from ions (cf. Naka- 
mura et al. 1997, Mahadevan & Quataert 1997), which effect partly 
reduces the differences between the high and low a models. These 
effects, and their impact on the radiative properties of the flow, are 
discussed in detail below. 

Figs[TJf) and[2tf) show also the dissipative heating rate of ions, 
2diss. The stabilized rotation results in a much stronger dissipative 
heating in the inner region for a = 0.998, however, this property is 
rather unimportant for effects investigated in this paper (neglecting 
both the direct heating of electrons and hadronic processes). 



3.1 Initial solutions 

We discuss here in some details the initial solutions although they 
are not self-consistent, because most of models available in litera- 
ture use the same local approximation for Comptonization. More- 
over, we amend some conclusions regarding the accuracy of local 
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Figure 1. Compaiison between initial and final solutions in models with M = 2 X 1 Mq and ;ji = 0. 1 . Panels (a-d) show the electron heating and cooling rates 
for the initial (a,c) and final (b,d) solutions. Q denotes a vertically integrated rate, so QR^ gives cooling/heating rates (per unit volume) times volume. The upper 
(a,b) and lower (c,d) panels are for a = 0.998 and a = 0, respectively. The short dashed curves show the synchrotron (magenta online) and bremssfi'ahlung 
(green) emissivities. The dot-dashed (black) curves show the Coulomb rate, Aie. The long dashed (red) curves show the effective heating rate of electrons, 
Aie — Gadv The dotted (cyan) curves show the compressive heating rate, Acompr- The solid (blue) curves show the global Compton cooling rate, Qcompt- The 
dotted (orange) curves in panels (a,c) show the Compton cooling rate in the local slab approximation, 2compt.loc ■ Panel (e) shows the ion (upper curves) and 
electron (lower curves) temperatures. The solid (black) and dashed (red) curves coirespond to our final, self-consistent solution for a = 0.998 and a = 0, 
respectively {T\ shown here is for both the initial and final solution); the dotted (black, a = 0.998) and dot-dashed (red, a = 0) curves show T^ for the initial 
solutions. Panel (f) shows the dissipative heating rate of ions (upper curves) and the Coulomb rate between ions and electrons (lower curves) for a = 0.998 
(solid, black curves) and = (dashed, red). 



approximations, derived previously in XIO. We find the initial so- 
lutions using the prescription for the Compton cooling rate in a slab 
geometry with initial photon energies of 1 eV, given in Dermer et 
al. (1991; hereafter D91). In XIO, we found that it gives a reason- 
ably good approximation for the innermost region. For the models 
considered here, the slab approximation appears to be less accurate, 
however, our choice of the slab case allows for a direct comparison 
with previous studies, a number of which used the slab approxi- 
mation based on D91, or its modification introduced in Esin et al. 
(1996). 

Figs[TJa,c) and|2la,c) show the heating and cooling rates in 
the initial solution. Comptonization of synchrotron photons is the 
most efficient radiative cooling process. However, the synchrotron 
emissivity decreases rapidly beyond its transition radius, r^ (^ 20 
for m = 0. 1 and ^ 100 for m = 0.5). Then, in the local model, 
the radiative cooling at c > r, is dominated by the (much weaker) 
bremsstrahlung and its Comptonization. Then, two regions can be 
distinguished in terms of the electron energy equation (ions are al- 
ways advection dominated), (i) At r > r^, the flow is adiabatically 
compressed, which is reflected in the increase of T^ with decreasing 
r (see Figs[T^ and|2^); both the radiative cooling and the Coulomb 



heating are much weaker than the advective terms, (ii) At r < r„ the 
radiative cooling becomes efficient and T^ decreases; the Coulomb 
heating is more efficient in this region and it exceeds the compres- 
sion work in the innermost part (except for the model with a = 
and m = 0.1). In this region, the advection of the internal energy 
contributes also to the heating of electrons, however, this effect is 
rather weak. 

The solid curves in Figs[Tfa,c) and[2ta,c) show the GR global 
Compton cooling rate obtained in our initial MC simulations, i.e. 
with the Te satisfying the energy equation ^ with the Compton 
cooling rate, !2compt,ioc, given by the local slab (D91) prescription. 
As we can see, there are significant deviations between the local 
(dotted curves) and global cooling rates, resulting from several ef- 
fects. Two of these effects should be generic to models of black- 
hole flows regardless of specific values of parameters of the flow. 
First, the input of seed photons from the inner region strongly en- 
hances the Compton cooling beyond r^. The second one involves 
the presence of the event horizon (neglected in the local approxi- 
mation). The presence of an inner boundary of the flow at the event 
horizon results in an obvious difi'erence with respect to the semi- 
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Figure 2. Same as in Fig. 1 but for M = IOMq and m = 0.5. 



infinite slab; namely, a large fraction of photons generated in the 
innermost region is captured, reducing the input of seed photons. 

Another effect affecting the structure of the flow is the depen- 
dence of the optical depth from a given point on direction. This 
strongly depends on the flow parameters and assumptions about 
heating and outflow. In particular, the flows for our parameters 
have the optical depth from a point in the equatorial plane along 
the outward radial direction, t^jju,, similar to the vertical one, T-, at 
any r, whereas r, „„, » t- in the model of XIO. Furthermore, for 
a = 0.998, both Tr,out and T; are much smaller than the optical depth 
in the inward radial direction, Trj„, except for the very innermost re- 
gion; in these models t^ju - T,,oi,t (- 0.2 for m = 0. 1 and == 1 for 
m = 0.5) at r » 2. For a = 0, T,-in == Tr,out (- 0.05 for ra = 0.1 
and - 0.3 for m = 0.5) at r » 10; in these models the difference 
between T,,i„ and T,-_out typically does not exceed a factor of 2. Then, 
for a non-rotating black hole, the local radiative properties would 
be more accurately described by the spherical geometry, while for 
a rapidly rotating black hole the local properties are intermediate 
between the sphere and the slab. Then, the assumption of the slab 
geometry overestimates the initial cooling rate in all models con- 
sidered here. 

Finally, as we note in XIO, the D91 formula is accurate for 
only two values of seed photon energies, 1 eV and 1 keV, and we 
use here the lower one. With this value, the formula gives a rea- 
sonable approximation for flows around a stellar black hole, where 
seed photons have energies > 1 eV. For flows around a supermas- 
sive black hole, however, with typical seed photon energies of 10"- 
eV (for M = 2x 10** Mq), it underestimates the cooling rate. Specif- 
ically, for this M we found deviations by a factor of ^ 3 between the 



formula of D91 for 1 eV and the results of our local slab MC sim- 
ulation using the exact energy distribution of seed photons (given 
by the synchrotron and bremsstrahlung emission distributions; see 
XIO for the description of our local MC model). 

3.2 Self-consistent solutions 

The initial solutions are qualitatively similar at both values of in 
considered here, the most significant dependence on m involves the 
increase of the value of rs with increasing ih. In contrary, detailed 
properties of the self-consistent, final solutions are more sensitive 
to ra. 

Figs[nb,d) show the heating and cooling rates of the final so- 
lutions for ra = 0. 1. At r > 10, a significant fraction of the power 
heating electrons is stored as their internal energy rather than ra- 
diated away. This is qualitatively similar to the initial solution, al- 
though 2rad is larger at r > 10 in the global cooling model and 
hence Qi„f is smaller than in the local cooling model; therefore, the 
temperature gradient is smaller in the outer region in our final so- 
lutions (see discussion below). We note that the strong contribution 
of the 2in( term significantly complicates the procedure for find- 
ing the self-consistent solution at this lower m, as both the global 
Compton cooling and the electron advection are strongly sensitive 
to, and affect, the electron temperature profile. 

Fig. [He) shows the electron temperature profiles for the ini- 
tial and the final solutions. We see that, in the final solution, T^ 
increases compared to the initial solution, by a factor of < 1.5 in 
the inner region, and up to a factor of 3 at several hundred ^g. The 
increase of T^ atr < 20 is an obvious effect, related with the over- 
estimation of the cooling rate by the slab approximation, combined 
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Figure 3. The radial profiles of hot flow parameters for M = 2 X 10^ Mq and m = 0.1. In all panels the solid (black online) and dashed (red) curves ai'e for 
a = 0.998 and a = 0, respectively, (a) The vertical optical depth, r, = ne(r,0)<TfH(nl2f-^. (b) The ratio of the scale height to the radius, (c) The azimuthal 
velocity, (d) The radial velocity. These parameters are determined in our initial (i.e., with the local slab approximation) solutions and are assumed to be the 
same in the final solutions. 




Figure 4. Same as in Fig. 3 but for M = 10 Mq and m = 0.5. 



with - crucial in that region - capturing of seed photons by the black 
hole, as discussed in Section [3TT] Somewhat surprisingly, the global 
solution requires higher T^ also at r > 20, where the global radia- 
tive cooling is significantly enhanced. The temperature must re- 
main high in this outer region because the energy balance could not 
be achieved for T^ decreasing too rapidly. Namely, we found that 
such a (too rapid) decrease of temperature at r > 20 would yield a 
large positive 2i„, term in the electron energy equation l[3]l, which, 
added to the (increased in global cooling) Q^j^a term, could not be 
balanced by the heating terms. We emphasise that this property is 
related to strong contribution of advective terms and it should not 
necessarily occur in the whole range of parameters, in particular at 
higher accretion rates. 

For m - 0.5 (Fig. |2j, the self-consistent T^ is higher in the 
innermost part than T^ in the initial solution, due to the same rea- 
sons as for m = 0.1. However, the Compton cooling is relatively 
more important in the outer region, which results in an approxi- 
mately isothermal structure at r > 4; in particular, we do not find a 
local maximum of T^, marking the transition (at — r^) from an adi- 
abatically compressed flow to an efficiently cooled flow in global 
solutions with m = 0. 1 and in all our local solutions. The temper- 
ature at which the energy balance is achieved is twice lower for 
a = 0.998, due to a much larger input of seed photons from the 
central region. 

The strong contribution of the compressive heating has cru- 
cial consequences for the radiative efficiency, ;; = L/Mc-, where 
L is the total observed luminosity of the flow, and for its scaling 
with m. We compare here our models with the same a for differ- 



ent m, although we note that they have also different M, which 
may slightly affect such a comparison (the main dependence on the 
value of M concerns the energy of synchrotron photons, influenc- 
ing the cooling rate). Comparing our results for ih = 0. 1 and 0.5, we 
find that the fraction of an accretion power transferred to electrons 
via Coulomb collisions varies approximately as Ai^^toi/Mc^ oc rii 
for a = 0.998 and for a = we note even stronger dependence on 
m. Then, heating dominated by the Coulomb coupling would lead 
to an approximately linear scaling of the radiative efficiency with 
m, as derived in previous works on this subject (see, e.g., eq. 12 
and references in Narayan & McClintock 2008; see also table 3 in 
Rajesh & Mukhopadhyay 2010). However, Aco„,p,-,tot/Afc^ is much 
less dependent, moreover negatively, on m; specifically, it is higher 
by a factor of ^ 3 for m = 0.1 than for m = 0.5. As a result, the 
dependence of 77 on m is significantly reduced. 



For a = 0, A, 



-compr.tot 



» Aie lot at m = 0. 1 and Ar, 



,>Aie 



at m = 0.5. The decrease of A, 



/Mc'^ with increasing m is par- 



compr.tot 

tially balanced by the increased contribution of Coulomb heating. 
Two further effects decrease the radiative efficiency at lower in: (1) 
the increase of the fraction of the energy given to electrons which is 
not radiated away (and increases the internal energy advected with 
the flow), and (2) the increase of the fraction of photons captured by 
the black hole (Section [333). As a result, the radiative efficiency 
for fl = is approximately the same for both considered values of 
111. More specifically, r] a 0.0042 (with L = 1.1 x 10"*^ erg s"' and 
L/Le = 4.2 X lO^") for m = 0.1 and 77 « 0.0045 (L = 2.9 x 10'"' • 



and L/Le = 2.3 x 10^^^) for m = 0.5. 
For a = 0.998, Aietot exceeds Acoi 



erg 



at both m = 0.1 and 
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Figure 5. (a) The angle-averaged spectra from the MC simulations for the 
initial (thinner curves) and final (thicker curves) solutions in models with 
M = 2 X 10** M0 and m = 0.1. The solid (black) and dashed (red) curves 
are for a = 0.998 and a = 0, respectively, (b) The angular dependence of 
spectra for the final solution with a = 0.998; the specfi'a from bottom to top 
are for/iob.s = 0.9-1, 0.5-0.6 and 0-0.1. 



0.5. Radiative efficiencies at these m differ only by a factor of 2, 
due to relatively stronger contribution of compressive heating at 
m = O.l and very strong reduction of the observed L by GR effects 
(Section [333J at m = 0.5. Specifically, for m = 0.1, ?/ ~ 0.0078 
(L = 2 X lO"*' erg s"' and L/Le = 7.8 x lO^") and for m = 0.5, 
77 a 0.016 (L = 1.1 X 10" erg s"' and L/Le = 8.3 x lO"'). We 
note that for a = 0.998 the apparent luminosity depends on S^bsi 
see Section l3.3.4l and the above values of L are averaged over all 
viewing angles. 

Figs |5ja) and [6ja) show the angle-averaged spectra for the 
initial and final solution. In all models, the normalization of the 
spectra of the final solution is much higher than that of the initial 
solution, due to the increase of the synchrotron emissivity result- 
ing from the increase of T^. The X-ray spectra are relatively hard, 
with the photon index F x 1.6, for both (in = 0.1, a = 0.998) and 
(m = 0.5, a = 0) and r K 1.5 for (ra = 0.5, a = 0.998); the last 
value is determined in the 2-10 keV range, as at higher energies 
a rather pronounced Wien peak leads to further hardening of the 
spectrum. For {m = 0.1, a = 0), the relativistic T^ and small opti- 
cal depth result in pronounced scattering bumps seen in this model 
spectrum; the average slope in the 2-200 keV range is F a 1.7. 
Thus, we see that the spectra harden with the increasing luminos- 
ity, for rii increasing from 0.1 to 0.5, from F ^ 1.7 to ^ 1.6 at a = 
and from F ^ 1.6 to ^ 1.5 at a = 0.998. 

In all models, the EL^ spectra have the maxima at several hun- 
dred keV. The electron temperature reaches the maximum values of 
~ 1 MeV for m = 0. 1 and ~ 400 keV for in = 0.5. 

Note, however, that although these flows produce photons with 



energies > m^c , we do not expect an efficient e* pair production 
in photon-photon collisions. The probability that a y-ray produces 
a pair is approximately given by the compactness parameter for 
this process, £{R) = Lx,RO'yy/(4mn^c^R), where Lx^r is the radiative 
power in photons with £ > 10 keV produced within R. It is « 1 
for all R in all models except for in = 0.5 and a = 0.998, where 
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Figure 6. Same as in Fig.|5]but for M = 10 Mq and m = 0.5. 



f > 1 in the innermost few /?„. The same conclusion was reached 
in previous studies of pair equilibria in advection-dominated flows 
(e.g. Esin 1999) where additional pair creation processes, namely 
electron-electron and photon-particle collisions, were also taken 
into account. 



3.3 Effects of black hole rotation 

The nature of the space-time metric (specified by the value of a) 
affects the radiative properties of hot ffows through various effects. 
As noted above, the stabilized circular motion results in a much 
larger density in the innermost region in models with large a. This 
affects the radiative efficiency in a manner discussed in Section [X2l 
here we restate some of these results with the emphasis on the de- 
pendence on a. Two further effects, discussed in Sections |3. 3. SI - 
IB. 3.41 concern the properties of photon motion in curved space- 
time. 



3.3.1 Heating 

Higher density implies a stronger Coulomb coupling for large a. 
For a = 0, Aie,tot is smaller than for a = 0.998 by a factor of 200 
and 35 for m = 0.1 and 0.5, respectively. The difference between 
the heating of electrons for high and low a is, however, strongly 
reduced by the high efficiency of compression, which shows very 
weak dependence on the value of a. Then, the total heating, Aje.tot + 
Acorapr,tot> 's higher in models with a = 0.998 only by a factor of 3.5 
(m = 0.1)and 16(m = 0.5). Extrapolating our results we can expect 
that at m <c 0. 1 the difference between the heating efficiencies in 
high and low a models vanishes. 

On the other hand, a much larger diff^erence between the heat- 
ing rates for large and small values of a could be expected in models 
involving the direct transfer of the dissipated energy to electrons; 
compare Aje and Acompi with 2diss shown in Figs[TJf) and[2jf). 
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3.3.2 Optical depth 

Higher density implies also a larger optical depth for large a. Then, 
despite a stronger heating in these models, the energy balance in 
the innermost region is achieved at smaller electron temperatures 
than for small a. On the hand, the highly efficient Comptonization 
occurs at very small r, therefore, a much more centrally concen- 
trated emissivity of Comptonized radiation is achieved in models 
with a = 0.998 and the observed luminosity is subject to stronger 
reduction by GR effects, discussed below. 



3.3.3 Collimation and capturing of radiation 

Different velocity fields result in a different strength of collimation 
toward the black hole horizon. In principle, larger values of the 
radial velocity, |v' |, imply a stronger collimation for a = 0, which 
could give rise to a larger fraction of photons being captured by 
the black hole. However, the radiation is produced, on average, at a 
smaller radial distance for a = 0.998 and, therefore, the reduction 
of the observed luminosity due to GR effects, including the photon 
capture and the gravitational redshift, turns out to be larger in that 
case. Specifically, in models with a = 0.998, Z^/Qrad.tot ~ 0.2 for 
in = 0.5 and 0.3 for m = 0. 1 ; the large reduction for m = 0.5 results 
from T > 1 (and trapping of photons) in the central region. For 
a = 0, L/2rad,iot ~ 0.7 (for in = 0.1) and 0.9 (for in = 0.5; in this 
case the emissivity profile is flatter). 



3.3.4 Intrinsic anisotropy 

For high values of a, effects unique for the Kerr metric, i.e. bend- 
ing of photon trajectories to the equatorial plane combined with a 
dependence of the gravitational shift of energy on the direction of 
photon escape, result in an intrinsic anisotropy of radiation pro- 
duced within the innermost few R„ (see Piran & Shaham 1977, 
Niedzwiecki 2005). The Comptonized component from that region 
is softer (i.e. has a larger spectral index and a smaller cut-off en- 
ergy) and has a smaller normalization when observed at a smaller 
Sobs. The gravitational effects, underlying this property, are strong 
only within 4Ag, then, the magnitude of the anisotropy in the total 
spectrum depends on the proportion between contributions from 
the innermost and the surrounding (beyond 'XRg) regions. 

Figs nib) and|6lb) show the t/^bs -dependence of the observed 
spectra in our models with a = 0.998 (in models with a = the 
dependence on 9^^,^ is negligible). As we can see, the spectra in- 
deed show some anisotropy but the effect is relatively moderate. 
For m = 0. 1, the total flux emitted face-on (//o(,s = 0.9-1) and edge- 
on (/Jobs = 0-0.1) corresponds to the isotropic luminosity of Ix 
and 3 X lO'*' erg s"', respectively. In this model, the compressive 
heating results in a rather flat emissivity profile, then, mixing of 
relatively weak contribution from r < 4 with radiation produced 
at more distant regions washes out the anisotropic properties. For 
m = 0.5, the face-on and edge-on fluxes correspond to 6.6 x 10^* 
erg s ' and 1.5 X 10^^ erg s' , respectively. In this model, the emis- 
sivity is strongly centrally concentrated; however, the optical depth 
of the inner region is large and a large fraction of photons produced 
there get reprocessed at larger distances or captured. At both ra, the 
spectral indices differ by AF ^ 0. 1 between a face-on and an edge- 
on observer. 

We expect that a much larger anisotropy would be produced in 
a model with moderate m (yielding r < 1) and strong direct heating 
of electrons (giving a very centrally concentrated emissivity). 



3.4 Illumination by an outer disc 

An optically thick disc may extend to relatively low radii in some 
objects observed at L/Le similar to these characterising our mod- 
els. The best example concerns XTE Jl 1 18+480, in which the de- 
rived truncation radius of the outer cold disc is - lOOiJg (Esin et 
al. 2001). We consider here the cooling of the inner hot flow due 
to Comptonization of thermal photons emitted by such a surround- 
ing, cold disc. For this we take the innermost part, within r = 100, 
of our self-consistent solutions described in Section 13.21 We use 
our MC method to compute the Comptonization of the seed pho- 
tons emitted by a cold, Keplerian disc extending from n,, = 100 to 
''out = 1000. The local black-body temperature of its thermal emis- 
sion is found using Page & Thorne (1974); note that this approach 
gives the largest possible luminosity of the outer disc as it involves 
the assumption that the rate of the outward transport of mechan- 
ical energy, which is the dominant source of the radiated energy 
at r > 100, is the same as in an untruncated Keplerian disc. We 
do not look for self-consistent solutions including this additional 
source of seed photons. Instead, we compute the total rate of the 
Compton cooling by the outer disc photons, Gcomnttof ^'^^ compare 
it with the total Compton cooling rate due to the synchrotron and 
bremsstrahlung emission, 2compt,iot! of our self-consistent solutions. 
We find that in models with a = the effect is small, with 
0.2, and for a = 0.998 it is negligible. 



iiCompt.tot' t^Compt.tot 



Q] 



JQc 



< 10 -. We note that the weakness of this ef- 



Compt.tol' i^Compt.tot 

fects results from the following properties. First, the total luminos- 
ity of the disc truncated at Hn = 100 is relatively small (specifically, 
it is approximately equal to the luminosity of the hot flow). Second, 
the solid angle subtended by the flow as seen by the outer disc is 
rather small. Third, the optical thickness of the outer parts of the 
flow is very small at the values of m considered here, so a small 
fraction of photons illuminating the flow get scattered. 

Our results are consistent with the conclusion of Esin et al. 
(2001) that the outer disc is not an important source for Compton 
cooling in this range of parameters. Obviously, the effect would 
be more important for a smaller truncation radius (cf. Esin 1997), 
at which the luminosity of the disc would be larger, as well as at 
higher m, at which a larger fraction of photons irradiating the flow 
would be scattered. 



3.5 Comparison with previous works 

Most of the specific properties discussed in this paper have been 
studied in previous works. However, their quantitative importance 
could not be assessed properly as neither of previous studies con- 
sidered self-consistently all relevant effects. 

The solution of a hot-flow structure based on a self-similar 
model (e.g., Narayan & Yi 1995) is commonly used due to its rel- 
ative simplicity. However, this simple solution introduces inaccu- 
racies of up to an order of magnitude even for the basic flow pa- 
rameters (e.g., the density given by eq. 1 in Mahadevan & Quataert 
(1997) is ~ 10 times larger than the density in our solutions), which 
leads to further inaccuracies in papers based on this model; e.g., we 
find that compressive heating dominates over Coulomb heating at 
much higher accretion rates than the critical value given by eq. (35) 
in Mahadevan & Quataert (1997). 

Gammie & Popham (1998) and Popham & Gammie (1998) 
presented an extensive study of various aspects of GR hydrodynam- 
ical description of a flow in the Kerr metric and they show that the 
black hole spin has a large effect on the velocity field, density and 
temperature of the flow (in ways also discussed in the present pa- 
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per); however, their study neglected radiative processes, so a more 
detailed impact on the observed spectrum could not be established. 

Manmoto, Mineshige & Kusunose (1997) included an ex- 
plicit description of radiative processes with local approximation of 
Compton cooling in a pseudo-relativistic model. The same model 
including outflows and direct electron heating was used by Yuan 
et al. (2003, 2005). Then, MOO improved such models by includ- 
ing a fully GR treatment of the hydrodynamical processes (but still 
with local approximation of Compton cooling). We have fully re- 
produced their results in our initial solutions; the accuracy of these 
models is directly illustrated by comparison of our initial and final 
(global) solutions. 

In an alternative (to MC) method for the treatment of global 
Comptonization in an optically thin flow, developed by Narayan, 
Barret & McClintock (1997), the flow is divided into a set of 
nested spherical shells and the iterative scattering method is used 
to describe their mutual interactions. The model was used, e.g. by 
Esin, McClintock & Narayan (1997), Quataert & Narayan (1999) 
and Esin et al. (2001). Neither of these papers presents any de- 
tails which could allow us to compare our cooling rates and only 
the produced spectra can be compared. The results of Esin et al. 
(2001), who use the GR hydrodynamical model from Gammie & 
Popham (1998) and include the special relativistic effects and the 
gravitational redshift in the calculation of radiative processes, is of 
particular interest for such a comparison. It appears that the model 
of Narayan etal. (1997) yields a higher cooling rate than our Monte 
Carlo model, as the spectrum computed in Esin et al. (2001) for pa- 
rameters similar to these in our models (see the solid curve in their 
fig. 1) has a break energy at ~ 100 keV, significantly lower than our 
model spectra. We do not know the cause of this diiference. 

Among previous studies of hot flows, the fully GR ADAF 
model of Kurpiewski & Jaroszynski ( 1 999, 2000), with a MC com- 
putation of global Comptonization, seems to be the closest to our 
present work. Their spectra show a much stronger dependence on 
a, see e.g. figure 4 in Kurpiewski & Jaroszynski (2000), which may 
be due to their values of m significantly smaller than in our models. 

Yuan et al. (2009) used an iteration method, similar to ours 
but with a diiferent approach to computing global Comptonization, 
to find the radiative cooling rate mutually consistent between the 
dynamics and radiation of the flow in a non GR model including 
strong outflows and electron heating. Taking into account signifi- 
cant differences between the results of our present work and XIO, 
pointed out below, we speculate that GR flows with a different 
strength of these two elfects may be characterised by different criti- 
cal accretion rates, and different maximum luminosities, than those 
found in the model of Yuan et al. (2009). 

Finally, we compare our results with the pseudo-relativistic 
model of XIO. First, we note that models based on the pseudo- 
Newtonian potential fail to properly describe the innermost region 
even for the a = case. Specifically, the pseudo-relativistic model 
predicts a decrease of the optical depth at r < 7 (see the dot-dashed 
curve in fig. 2b in XIO) while the fully GR model with the same 
parameters (see our Fig. [3^) yields the optical depth monotonically 
increasing with decreasing r. The reason for this discrepancy can 
be traced to the (unphysically) large radial velocity in the pseudo- 
relativistic model, formally exceeding c at small r, which implies a 
much smaller density than that in the GR model. 

Second, in the present work we use the equation for vertical 
equilibrium in the form derived for the Kerr metric by Abramowicz, 
Lanza & Percival (1997). In XIO we used its simplified version, 
H = c's/flK, where c, is the sound speed and Q.k is the Keplerian 
angular velocity, and we note significant differences in scale heights 



derived in these two approaches. Specifically, the equation from 
Abramowicz et al. (1997) yields a larger scale height, up to a factor 
of ~ 2, for the flow parameters describing our models. 

Third, additional effects included in XIO, namely the strong 
outflow and direct viscous heating of electrons affect the flow struc- 
ture in a complex manner and lead to a qualitatively different effect 
of taking into account the global nature of the Comptonization pro- 
cess. Obviously, the outflow leads to a reduction of t- at small r. 
Therefore, Tj is much smaller, by a factor of ~ 10 at r < 10, in the 
model of XIO than in our present model with a = and m = 0.5 
(which parameters are the same as those in XIO). Then, the outflow 
reduces the compressive heating of ions via reducing the density 
gradient, and the direct heating of the electrons reduces the viscous 
heating of ions. Also, the outflow reduces the viscous heating rate 
per unit volume, flattening its radial profile, but the viscous heat- 
ing rate per ion remains unchanged. We note that the compressive 
heating of ions exceeds their viscous heating in our models, and it 
is much larger, by a factor of TJT^, than the compressive heating 
of electrons. Together, both reductions of the ion heating cause a 
reduction of Tj, which, in turn, leads to a decrease of H/R by ~ 30 
per cent (with respect to the pseudo-relativistic model without both 
an outflow and electron heating) within the innermost several R„. 

The simplified hydrostatic equilibrium condition together with 
the outflow and electron heating considered in XIO yield a much 
smaller H/R than our present model (compare our Fig.|4j3 with fig. 
3b in XIO). As a result, in spite of the strong mass loss to the out- 
flow, the total radial optical depth in the model of XIO is still higher 
than that in our model with a = and m = 0.5. Also, the radial op- 
tical depth of XIO is much larger than t.; then, the geometry of the 
flow resembles that of a slab. 

Fourth, comparing the radiative efficiencies of the XIO model 
and the present one with a = 0, m = 0.5 illustrates the strong impact 
of the direct heating of electrons. Namely, in XIO the value of r] is 
three times larger, although most of the potential energy of the flow 
is lost to the outflow in their model. 



3.6 Strength of magnetic field and electron heating 

We discuss here briefly our assumption on the strength of the mag- 
netic field as well as our neglect of the viscous heating of elec- 
trons. Our assumption of rather weak magnetic field, with the mag- 
netic pressure of 0. Ip, is supported by results of the magnetohydro- 
dynamic (MHD) simulations in which amplification of magnetic 
fields by the MRI typically saturates at such a ratio of the mag- 
netic to the total pressure (e.g. Machida, Nakamura & Matsumoto 
2004 and references therein). It is rather easy to assess changes 
of the flow structure and the spectrum resulting from a change of 
this assumption. Namely, a stronger magnetic field would result in 
a stronger synchrotron emission and this would cause a decrease 
of the electron temperature and a softening of the Comptonization 
spectrum (see, e.g., Esin et al. 1997). 

An issue of the direct viscous heating of electrons is more un- 
certain. Current applications of ADAF models typically assume 
large values of electron heating to total heating, 6 ~ 0.5 (e.g.. 
Yuan et al. 2003), however, the model has a degeneracy between 
the value of S and the strength of an outflow (Quataert & Narayan 
1999). Various attempts to assess this effect on theoretical grounds 
seem to favour small values of (5 for weak magnetic fields and rather 
high in, which are considered in our models. Analytic investigation 
of particle heating by MHD turbulence in ADAFs indicates that 
the turbulence primarily heats protons for weak magnetic fields, 
while electrons are primarily heated for strong magnetic fields (see 
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Quataert & Gruzinov 1999). However, large uncertainties in the di- 
vision between electron and proton heating result from uncertain- 
ties in the description of the turbulence; ^63 = 0.9, assumed in our 
models, may correspond to the values of (5 between =^ 0.006 and 0.6 
(see fig. 2 in Quataert & Gruzinov 1999). Details of another poten- 
tial mechanism of electron heating discussed by Bisnovatyi-Kogan 
& Lovelace (1997) and Quataert & Gruzinov (1999), i.e. reconnec- 
tion, remain even more uncertain. 

Simulations of MRI turbulence by Sharma et al. (2007) in- 
dicate that pressure anisotropy, created in turbulent plasmas, may 
give an additional mechanism for particle heating; they approxi- 
mate the fraction of the viscous energy which heats electrons as 
S ~ {T^/Ti)'-'^/3. However, their model assumes a fully collision- 
less plasma, which approximation is valid only for m « 0.06. At 
higher accretion rates. Coulomb collisions suppress the electron 
pressure anisotropy and electron heating is negligible. Then, their 
prescription for 6 seems to be not relevant for our models. 

We note that the direct viscous heating would exceed the 
Coulomb and compressive heating of electrons for S^O.l in mod- 
els with a = 0, and for 5 > 0.01 in models with a = 0.998. For 
these values of S, flows should have larger luminosities than these 
assessed in our models. We can also expect that flows with such val- 
ues of S produce harder X-ray spectra, with higher cut-off energies, 
than obtained in our models. However, we note also that large val- 
ues of S, for which the heating of ions is significantly reduced, may 
have a strong impact on the flow structure, especially if a strong 
outflow is also present, which may lead to enhanced global Comp- 
ton cooling (as discussed in Section [331 ). Then, the above simple 
estimation of spectral changes related to the increase of 6 may be 
incorrect; in particular, we note that the spectrum obtained in our 
model with S = 0.5 in XIO has a smaller cut-off energy than our 
spectra computed here with 6 = 0. 



4 COMPARISON WITH OBSERVATIONS 

4.1 FR I radio galaxies 

A substantial observational evidence indicates that accretion flows 
typically proceed through a radiatively-inefficient mode below a 
characteristic luminosity of == O.OlLg in AGNs (e.g.. Ho 2008; see 
Yuan 2007 for a review on the applications of hot flow models to 
low luminosity AGNs). FR I radio galaxies make an interesting 
class of objects observed at such luminosities. They do not show 
signs of the presence of an optically thick accretion disc (e.g., Chi- 
aberge, Capetti & Celotti 1999). Remarkably, an average radiative 
efficiency, rj < 0.005, estimated for FR Is by Balmaverde, Baldi & 
Capetti (2008) is consistent with the values of ;; obtained in most 
of our solutions. 

According to our results, hot flows with relevant accretion 
rates should produce rather hard spectra, extending to MeV ener- 
gies. Then, MeV observations are crucial for verification whether 
the high energy component is produced by the thermal Comp- 
tonization in a hot flow; however, the present data quality in this 
energy range is too poor for most objects to allow for such analy- 
sis. 

Thanks to its proximity, Centaurus A is the best studied FR I 
radio galaxy, suitable for detailed testing of various accretion mod- 
els through modelling of its broadband spectrum. Its X-ray lumi- 
nosity (Lo-io/Le ~ 6 X 10"', see below, where L2-10 is the 2-10 
keV luminosity) satisfies the criterion, suggested in Wu, Yuan & 
Cao (2007), for a FR I galaxy to be dominated by a hot flow rather 



than a jet emission, of L2-10 ^ 5 x IO^^Le. Indeed, many works at- 
tributed its spectral components to emission from accretion flow 
(e.g. Evans et al. 2004; Whysong & Antonucci 2004; Meisenheimer 
et al. 2007) although they can also be explained in terms of jet mod- 
els (e.g. Abdo et al. 2010). Below we briefly compare predictions 
of our model with the observed properties of the X-ray emission in 
CenA. 

(i) Luminosity. Over the last decade, Cen A was typically observed 
with the unabsorbed X-ray luminosity of L2-10 - 5 x lO"*' erg s"' 
(e.g. Evans et al. 2004, Markowitz et al. 2007). Our model with 
lit = 0.1 and a = 0.998 gives L2-10 = 1.4 x lO''^ erg s"' (assuming 
(?o(,s = 50°), which is approximately consistent with the average 
observed L2-10 taking into account recent estimates of its black- 
hole mass of ^ 6 x 10' Mq (e.g. Haring-Neumayer et al. 2006), 
which is ~ 3 times smaller than M assumed in our model. Also for 
lit = 0. 1 and a = we get a roughly consistent L2-10 = W*^ erg s"' . 

(ii) Spectral index. The time-averaged spectral index of the X-ray 
emission is F « 1.7-1.8 (Beckmann et al. 2011, Rothschild et al. 
2011). The model spectrum for a = 0.998 is slightly harder, F « 
1.6. The average slope for a = is roughly consistent with the 
observed one, but this model spectrum has pronounced scattering 
bumps, which have not been observed in Cen A. 

(iii) High energy cut-off. Cen A is the only FR I galaxy with data 
at MeV energies. Steinle et al. (1998) have found a spectral break 
in the F ^ 1.7-1.8 spectrum at ~150 keV Another power law spec- 
tral component was found above the break up to the MeV range. 
Grandi et al. (2003) and Beckmann et al. (201 1) fitted the hard X- 
ray spectra by an e-folded power law and found e-folding energies 
of ~500 keV, which is roughly compatible with the break energy of 
~150 keV. Rothschild et al. measured the power law spectrum up 
to 200 keV, within which energy range there was no indication for 
a break or cut-off. 

(iv) Variability. The spectral index remains roughly unchanged 
while L2-10 varies by a factor of 3 (Rothschild et al. 201 1). In our 
model, changes of L2-10 can be attributed to changes of a few differ- 
ent parameters, namely, the strength of the magnetic field, the frac- 
tion of accretion power directly heating electrons or the accretion 
rate. The results of our current study allow us to assess only effects 
related to changing of m (which case is most commonly considered 
as the explanation of luminosity variations), and we note that the 
model predicts a small change of the photon index, by AF =: 0.1, 
corresponding to the change of L/Lg by a factor 10. This appears to 
be compatible with the reported spectral variability within F ^ 1.7- 
1.8 only. 

In summary, our model can reproduce most of the properties 
of the X-ray emission in Cen A. However, it is inconsistent with 
the relatively low observed value of the break energy of ~ 200 keV 
Furthermore, the a = 0.998 model predicts a slightly harder spec- 
trum than that observed. However, an additional process, which 
would lead to a softening of the spectrum and a decrease of the 
cut-off energy, has to be included in the model with a - 1. 
Namely, an efficient nonthermal synchrotron emission from rela- 
tivistic e*, coming from the decay of charged pions copiously pro- 
duced through proton-proton collisions in a flow around a rapidly 
rotating black hole, where Ti > 10'^ K, would provide a strong seed 
photon input in addition to the thermal synchrotron emission. How- 
ever, we leave a more detailed comparison of the spectrum ob- 
served from Cen A with the hot flow model spectra to a future work, 
in which we will take into account the relevant hadronic processes. 
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4.2 Hard spectral states 

Black-hole binaries in the hard spectral state have typical X-ray 
photon spectra (measured typically in the 3-10 keV range) of F ^ 
1.5-2 (e.g., Zdziarski & Gierlihski 2004; McClintock & Remillard 
2006). In luminous spectral states, L>0.01Le or so, the high en- 
ergy cut-offs in those spectra occur at energies ^ 50-200 keV (e.g.. 
Grove et al. 1998; Wardzihski et al. 2002; Zdziarski & Gierliriski 
2004), as measured by the position of the spectral turnover in EFg 
spectra. 

Thus, the X-ray slope of our model spectra, T ^ 1.5-1.7, is 
in the range of the observed values. On the other hand, the high- 
energy turnovers in our spectra are at =^ 500 keV, which is higher 
than the values observed. However, our models have L/Le - (0.4- 
8) X 10"^, which is lower than L/Le at which most of the spec- 
tral cut-offs have been measured. When measured, the cut-off en- 
ergy has been found to increase with the decreasing luminosity 
(Wardziriski et al. 2002; Yamaoka et al. 2006; Miyakawa et al. 
2008), and thus it is possible that the cut-off energy is at ~500 keV 
in some low-L hard states. 

The value of the X-ray photon index shows a negative corre- 
lation with the luminosity in luminous hard states of black-hole bi- 
naries, at L > O.OILe or so, i.e., the spectra soften with the increas- 
ing L. However, the correlation changes its sign at low L, where 
the spectra soften with the decreasing L (e.g.. Corbel et al. 2004; 
Zdziarski et al. 2004; Yuan et al. 2007; Wu & Gu 2008; Sobolewska 
etal. 2011). 

As pointed out by Sobolewska et al. (2011), this indicates a 
change in the physical nature of the X-ray source at a luminosity 
of a fraction of a per cent. Those authors discuss two scenarios 
explaining this change. One is that the emission is dominated by 
an accretion flow at high L, and by a nonthermal jet radiation (most 
likely synchrotron) at low L. This scenario for the black-hole binary 
XTE J1550-564 is discussed by Russell et al. (2010), who found 
the near-infrared/optical (~ lO'"*^ Hz) emission is tightly correlated 
with X-rays at low L, but not at high L. A potential problem for this 
model is presented by the fact that basically the same L-T correla- 
tion, negative at low L/Le and positive at high L/Le is observed in 
AGNs. The transition is at L/Le ~ 0.01 or so, about the same as in 
black-hole binaries (e.g., Gu & Cao 2009; Constantin et al. 2009; 
Younes et al. 201 1; Veledina et al. 201 1). The magnetic field in jets 
in black-hole sources is expected to scale as 6 oc M"''^ (e.g., Heinz 
& Sunyaev 2003), which (since the cyclotron frequency is oc B) 
very strongly affects the jet emission, leading, in particular, to the 
strong mass dependence of the relative strength of the jet emission 
at a given frequency (Merloni, Heinz & Di Matteo 2003). Thus, 
it appears highly unlikely that the X-ray jet emission would start 
to dominate the accretion-flow emission below about the same Ed- 
dington ratio in both black-hole binaries and AGNs (see also Yuan 
& Cui 2005; Yuan, Yu & Ho 2009). 

The second possibility explaining the form of the L-T correla- 
tion discussed by Sobolewska et al. (2011), and independently by 
Veledina et al. (2011), is that it is due to a change of the charac- 
ter of the dominant seed photons for Comptonization. The overall 
idea follows the model for spectral changes of black-hole binaries 
developed by Esin et al. (1997, 1998), whose model did predict 
the changing L-T correlation, see, e.g, fig. 1 in Esin et al. (1998). 
At high L, above the transition L/Le, the seed photons are likely 
to be from an optically-thick disc overlapping with the hot flow. 
There is a lot of observational evidence for the presence of an 
outer cold disc in luminous hard states (as summarized by Done, 
Gierliiiski & Kubota 2007), e.g., a strong correlation of T with 



the relative strength of Compton reflection (Zdziarski, Lubiriski & 
Smith 1999). Also, the inner disc radius decreases with the increas- 
ing L (e.g.. Done et al. 2007). This naturally explains the positive 
L-F correlation, since the resulting increasing overlap between the 
cold disc and the hot flow leads to both an increase of the strength 
of Compton reflection and increased cooling of the Comptonizing 
plasma, resulting in softening of the X-ray spectra. However, at low 
L/Le, the outer disc may be truncated far away, and the dominant 
source of seed photons for Comptonization becomes synchrotron 
emission, as in the models studied in this work as well as in the 
ADAF model in general (Narayan & Yi 1995). This model can also 
explain the correlation of the X-rays with the near-infrared/optical 
emission, which, in this scenario, would be from the ADAF flow, 
see Fig.|6] 

As shown in Section [T2l our models spectra do harden with 
the increasing L, with AF ~ -0.1 for L increasing by a factor of 
several. This strongly supports the suggestion of Sobolewska et al. 
(201 1) that the negative L-F correlation seen at L/Le "Z ^-^^ t^ due 
to the dominance of synchrotron seed photons in that regime. A 
hardening of the ADAF spectra with increasing L, below a few per 
cent of Le, was studied previously by Esin et al. (1997); their model 
seems to predict a much stronger effect, with the increase of m by 
a factor of 10 corresponding to AF a -0.7. In this work, we have 
considered only two values of the accretion rate, and have not fully 
tested this model against observational data. We will consider this 
issue in detail in a paper in preparation. 

We note that Oda et al. (2010) propose another possible expla- 
nation for the difference between the low and high luminosity hard 
states. It is a transition between two equilibrium states of an opti- 
cally thin, two temperature flow with a weak and a strong magnetic 
field, respectively. Such a change of magnetic field is outside the 
scope of our model. 



5 CONCLUSIONS 

We have developed a fully GR model with global Compton scat- 
tering of hot, tenuous accretion flows. The dominant source of seed 
photons in our model is the synchrotron process. We assume that an 
outer, optically-thick, accretion disc is truncated at large radii, so ir- 
radiation of the hot flow by it is negligible. We have also neglected 
direct viscous electron heating and outflows. We have considered 
the accretion rates of m = 0.1 and 0.5, for which our model gives 
the bolometric luminosities of L/Le - (0.4-8) x 10"'. 

In this range of L/Le, we find our models predict X-ray spec- 
tra hardening with increasing L, which is in agreement with the 
correlation observed at L/Le 'Z 0-01 in both black-hole binaries in 
the hard state and AGNs. This is also in agreement with the sug- 
gestion of Sobolewska et al. (2011) that synchrotron seed photons 
dominate in this range of L/Le, whereas an optically-thick accre- 
tion disc overlapping with the hot flow provide dominant source of 
seed photons at higher L. 

We have compared our models to the low-L FR I radio galaxy 
Cen A. We have found a good agreement with the observed X- 
ray spectral slope, but the position of the high-energy break, ob- 
served at ~200 keV, is lower than the model cut-off energy. This 
may be possibly explained by the hadronic processes, which would 
be present if the black hole in this source has a high spin. 

We have studied the impact of the black hole spin and we find 
that despite a very strong influence of this parameter on the dy- 
namical properties of the inner part of the flow, the difference be- 
tween the observed radiation from a flow surrounding a maximally 
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rotating and a non-rotating black hole, with other parameters un- 
changed, is rather moderate. In particular, their luminosities differ 
only by a factor of 2-3 (i.e. smaller than a factor of == 6 in Keplerian 
discs). We find this is due to a strong role of the compressive elec- 
tron heating, which process only weakly depends on the black-hole 
spin. We note, however, that a much larger dependence on the spin 
value is expected in models involving direct heating of electrons, 
or taking into account hadronic processes. 

We stress that a proper treatment of global Compton scatter- 
ing is crucial for accurate modelling of the flow structure, most 
importantly, for the electron temperature profile, and for the spec- 
tral formation in optically thin flows. In models with m = 0.5, we 
find the entire inner flow is efficiently cooled by Comptoirization 
when the non-local nature of this process is taken into account. At 
m = 0.1, global Comptonization is important, but also electron ad- 
vection effects are strong at f > 10, giving rise to large temperature 
gradients. As compared to the results with the global Comptoniza- 
tion, the local model assuming a slab geometry provides in general 
a very poor approximation, overestimating the cooling rate by an 
order of magnitude in the innermost region. 
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